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Modern density functional theory (DFT) calculations employ the Kohn-Sham (KS) system of 
non-interacting electrons as a reference, with all complications buried in the exchange-correlation 
energy (-Exc). The adiabatic connection formula gives an exact expression for E xc - We consider 
DFT calculations that instead employ a reference of strictly-correlated electrons. We define a 
"decorrelation energy" that relates this reference to the real system, and derive the corresponding 
adiabatic connection formula. We illustrate this theory in three situations, namely the uniform 
electron gas, Hooke's atom, and the stretched hydrogen molecule. The adiabatic connection for 
strictly-correlated electrons provides an alternative perspective for understanding density functional 
theory and constructing approximate functionals. 



I. INTRODUCTION 

For most modern calculations using density functional 
theory (DFT) pQ, the accuracy of results depends only on 
approximations to the exchange-correlation functional, 
Eyic[n]. An exact expression for i? xc [n] is given by 
the adiabatic connection formula [2J |3], in which E^ c 
is expressed as an integral over a coupling constant A, 
which connects the reference (Kohn-Sham system, A = 0) 
and the real physical interacting system (ground state, 
A = 1), keeping the density n(r) fixed. Study of the 
adiabatic connection integral has proven very useful for 
understanding approximate (hybrid) functionals |3J [5], 
and is an ongoing area of research [HI El IE]- 

Almost all modern DFT calculations employ the Kohn- 
Sham (KS) system [3] as a reference. The KS system is 
defined as the unique fictitious system that has the same 
density as the real system, but consists of non-interacting 
electrons. The great practicality of KS DFT is due to the 
relative ease with which the non-interacting equations 
can be solved, with relatively crude approximations, giv- 
ing KS DFT a useful balance between speed and accu- 
racy. 

However, DFT calculations could also be based on an- 
other fictitious system which is known as the strictly- 
correlated (SC) system [TO]. The strictly-correlated sys- 
tem has the same density as the real system (as does 
the KS system), but the Hamiltonian consists of elec- 
tron repulsion and external potential energy terms only. 
In recent years, the pioneering work of Seidl and oth- 
ers QJU HU EE2 CE3 OS H3 US E] has led to substan- 
tial progress in solving this problem exactly and effi- 
ciently. The strictly-correlated electrons (SCE) ansatz 
[HI ITFJl [Pol ITT] has been shown to yield the density and 
energy of this system, going beyond the earlier point- 
charge-plus-continuum (PC) model [12] [13]. They have 
achieved great success in calculating spherical symmetric 
systems with arbitrary number of electrons |15j . 

In this article, we look to the future and assume that 
the strictly-correlated limit of any system can be cal- 
culated with less difficulty than the original interacting 
problem, and all our successive work is developed based 
on this reference. We derive a new version of the adia- 



batic connection formalism, which connects the strictly- 
correlated system (fully interacting) and the physical sys- 
tem. We also introduce a new coupling constant /i, and 
a "decorrelation energy" _E DC , the counterpart of E^ c in 
KS DFT, which must be evaluated to extract the true 
ground-state energy from the calculation of the strictly- 
correlated system. We argue that, as long as the strictly- 
correlated system can be solved easily (just as the KS 
case) , one can develop another version of DFT based on 
this system, a version that is better-suited to strongly 
localized electrons. 

Throughout this paper, we use atomic units (e 2 = ft — 
fx = 1), which means that if not particularly mentioned, 
all energies arc in Hartrees, and all lengths are in Bohr 
radii, etc. 



II. THEORY 

In this section, we introduce the alternative adiabatic 
connection formula, and relate its quantities to more fa- 
miliar ones. All results here are formally exact. 



A. Kohn-Sham Adiabatic Connection 

In KS DFT, the total energy for the interacting 
ground-state is expressed as: 

E[n] = T s [n] + J d 3 r t> cxt (r) »(r) + U[n] + E xc [n]. (1) 

In this equation, T s is non-interacting kinetic energy 
of KS orbitals that are eigenfunctions of the non- 
interacting KS equation, w cx t(r) is external potential (nu- 
clear attraction in the case of atoms and molecules), U 
is the Hartree energy defined as the "classical" Coulomb 
repulsion between two electron clouds, and i? xc is the 
exchange-correlation energy 18J. The adiabatic connec- 
tion integral pO[S] then gives an exact expression for E^ c : 

E xc [n}= f d\W x [n], (2) 
Jo 
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where one can show that W\ =< ^E' A |Ke|^' A > —U, in 
which ^ A is the wavefunction that minimizes T + AV^c 
but has the same density as the real ground-state system 
[TT?] . At A = 0, one recovers the KS system, and at A = 1, 
one recovers the real interacting system. In this way, one 
connects the KS system with the real interacting system 
by changing A from to 1. 

A cartoon of W\ versus A is shown in the upper panel 
of Fig. [I] By definition, we have Wo — E x and the area 
under the curve is E xc . We can also identify the kinetic 
correlation energy T c — E xc — W\ in this graph. 
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electron-electron repulsion energy alone, subject to re- 
producing the exact density. In practice, there might be 
multiple degeneracies, so it is best defined in the limit as 
the kinetic energy goes to zero. 

Then, using the strictly-correlated (SC) system as the 
reference, the energy of the true interacting ground state 
is: 

E[n) = U sc [n] + J d 3 rv cxt (r)n(r)+T s [n} + E DC [n], (3) 

where U sc =< I- 00 |t> cc |1' 00 >. In KS DFT quantities, 
U sc = U + W QO [see Eq. (§]. 

Just as we separate the Hartree energy from the total 
energy in KS DFT [Eq. 0], here in Eq. p]) we sepa- 
rate T s [n] from the total energy in SC DFT. There are 
a variety of algorithms that can be used to extract this 
quantity accurately for any given density, effectively by 
inverting the KS equations [20]. We label the remainder 
as the "decorrelation energy" , E DC [n] . The reason we call 
it "decorrelation energy" is that, if we consider the elec- 
trons in the reference system "strictly correlated", with 
energy [/ sc [n], the electrons in the real system are less 
correlated than in the reference system. We will see the 
physical meaning explicitly very soon. 

So far, we have defined our reference, and next we de- 
duce an exact expression for the newly-defined "decor- 
relation energy" i? DC [n] with the adiabatic connection 
formalism, just as one does for E xc [n] [TUl [H] in the KS 
DFT [Eq.(§]. 




FIG. 1: Traditional (upper panel) and the new (lower panel) 
adiabatic connection curves. 



B. Strictly-Correlated Reference System 

The KS wavefunction can be defined as the wavefunc- 
tion that minimizes the kinetic energy alone, but whose 
density matches that of the interacting system. Analo- 
gously, the SC wavefunction is found by minimizing the 



C. Strictly-Correlated Adiabatic Connection 
Formula 

We denote ^ as the wavefunction minimizing = 
H 2 T +V CC + V^t w ith density n(r). For fi = 0, we recover 
the strictly-correlated system, and for // = 1, we recover 
the real system. For each value of fi, there is a corre- 
sponding unique external potential yielding the correct 
density, 7^ t (r). So we have: 



V r , 



Using Hcllmann-Feynman theorem [21 , we have: 



(4) 



dE^ 





dH» 




dfj, 



2/if 1 ^ 



dfi 



(5) 

Integrating and cancelling the external potential terms at 
both sides, we recognize the left hand side is just T 3 [n] + 
E DC [n]. Thus: 



E DC [n] 



dfl2fi (\fr 



T 



T s [n]. 



(6) 



This is our adiabatic connection formula for strictly- 
correlated electrons. Finally, since T c [n] = T[n] — T s [n], 
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and T s [n] is independent of /i: 

E U c[n] = / d/iK^n], 
Jo 



where 



K„ 



2 M T£[n] 



(7) 



(8) 



This is the SC doppelganger of Eq. We plot a car- 
toon of the integrand vs. /i in the lower panel of 
Fig. [I] identifying the area below the curve as E DC , and 
noting that K x = 2T C . 



D. Relation to Kohn-Sham DFT 

From a formal viewpoint, what we derived here is not 
new, but simply another way to describe the real inter- 
acting system. Thus we can relate all quantities defined 
here, such as -B DC [n] and -ft^n], to quantities defined in 
the traditional KS DFT. Since both Eq. ^ and Eq. 
are exact for the real system, and if we use the expression 
of U sc [n] in KS language [see discussion below Eq. jjjj], 
we find: 



-EdcN = E xc [n] — Woo[n]. 



(9) 



Thus -B DC [n] defined in our theory is just the difference 
between the usual exchange-correlation energy of the real 
system, E xc [n], and the potential contribution to the 
exchange-correlation energy of the strictly-correlated sys- 
tem, Woo[n]. 

We can also deduce an expression for in terms of 
W\. Since ^ minimizes = fi 2 T + V cc while yielding 
n(r), and ^> x minimizes T+XV CCl we deduce vff 1 /^ 
Now, from the scaling properties of KS DFT [! 
know: 



, we 



dE x 



(10) 



If we write = T> 
tribtion to £ C A , = X{W, 



dTl 
d\ 



x 
X 



e., U£ is the potential con- 
fix): we have: 



(11) 



Integrating over A from to 1 //i 2 , and using the defini- 
tion of W\ in Eq. (|2| and that E x = A£" x by scaling [22] , 
we can express Tcjn] in terms of W^Jn], we find: 



2/z 



i/V 2 



dX(W x [n\-W x/lJ ?[n\)- (12) 



From this relation, we can generate the new adiabatic 
connection curve, as long as we know the integrand of 
the KS adiabatic connection, i.e. Wa[ti] for A = 1 to oo. 



E. Exact conditions 

Many of the well-established exact conditions on the 
correlation energy can be translated and applied to the 
decorrelation energy. In particular, the simple relations 
between scaling the density and altering the coupling con- 
stant all apply, i.e., 

E x [n] = X 2 E c [n x/x ], (13) 

where ?i 7 (r) = 7 3 71(71*). Thus, in terms of scaling: 

K^n] = \T c [n^]. (14) 

Note that, as ju — > 00, — > 0, while K^ =0 = IW'^, 
where W' is defined in the expansion of W x as A — > 00 



= lim 

A — >oo 



^(VFa-IFoo). 



(15) 



Thus the SC energy is found from solving the strictly- 
correlated system, while K^ = q is determined by the zero- 
point oscillations around that solution. Both are cur- 
rently calculable for spherical systems [TH1 HZ] . 

The most general property known about the correla- 
tion energy |22j is that, under scaling toward lower densi- 
ties, it must become more negative. In turn, this implies 
that W\ is never positive. Using the definition of T£ and 
changing variable A = 1/n 2 in Eq. (Ill, we find: 



dTg 
dfi 



2 dW\ 
/i 5 dX 



<0, 



(16) 



then using = 2/iTc and the fact that > 0, we 
find: 



dfi 



(17) 



Also, because T£ > 0, so = 2^iT§ > 0, and E DC , as 
an integration of K^, is always positive. 

Based on these properties of K^, a crude approxima- 
tion to can be a simple exponential parametrization, 
using Kq and the derivative of In K at fi — as inputs: 



7 = 



d[i 



■In if 



(18) 



III. ILLUSTRATIONS 



In this section, we illustrate the theory developed 
above on three different systems, to show how be- 
haves for very different systems, and where the adia- 
batic connection formula might be most usefully approx- 
imated. 
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A. Uniform Electron Gas 

For a uniform electron gas, we assume we know the cor- 
relation energy per particle, e c , accurately as a function 
of r s — (3/47rn) 1 / 3 . In order to apply Eq. (12 1 to calcu- 
late K u \n\, we use e^.(r s ) = A 2 e c (Ar s ) [55]. Substituting 
into Eq. (10), changing variables A = l//i 2 , and using 
= 2^rf~= Nk^, with N the number of particles, we 
find: 



„unif 



■> - ^3 drs ( r Mr s ))\ rs/fl 2 ■ 
Using Eq. ^ and the definition of W^, we find: 

^unif _ c 

£ DC " £ C1 ! 

where do is defined below and do = 0.433521. 
large- r s limit or the low-density expansion |23j : 



(19) 

(20) 
In the 



Cc^s) = h 



di_ 

3/2 



— 



-3.66151 from data of Ref. 123] 



where di 

ing this expansion into Eq. (19 1, we find: 
d 2 

r. 



(21) 
Substitut- 



K 



unif 



dl_ 

3/2 



as \i -> 0. (22) 



Thus /Cp is expected to have a well-behaved expansion in 
powers of \x for small /i. 

Using Perdew and Wang's [23 parametrization of the 
correlation energy of the uniform gas, we plot vs. fi 
for r s = 1 in Fig. [2] and find e DC = 0.374 at r s = 1. 



and 7 = 5.0826. We plot the exponential parametriza- 
tion in Fig. [2] and we can see that it decays much faster 
than the exact curve, producing a e DC that is too small 
by about 25%, which means about 150% larger in |e c | 
[see Eq. Q]. 

We calculated e D c/kc| for different values of r s and 
plot the curve in Fig. pi At small r s , e DC 3> |e c |, which 
suggests that the KS reference system is a better starting 
point, as a smaller contribution to the energy needs to be 
approximated. At large r s , |e c | > e DC so e DC is a smaller 
quantity and may be better approximated. Under such 
circumstances, the strictly-correlated system might serve 
as a better reference. For the uniform gas, the switch- 
over occurs at about r s = 16, which is at densities much 
lower than those relevant to typical processes of chem- 
ical and physical interest. However, as we show below, 
for systems with static correlation, this regime can occur 
much more easily. 
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tron gas (r s = 1) and a simple exponential parametrization. 



Using the exact curve for r s = 1 in the simple expo- 



FIG. 3: £dc/|£c| for different r a for uniform electron gas. 



B. Hooke's Atom 

As we pointed out, as long as we have an approximate 
WaW for A between 1 and 00, we can substitute it into 
Eq. ( 12 ) to get the new adiabatic connection formula for 



nential parametrization [Eq. ( IS l] , we find kq = 1.44073 



the decorrelation energy. Of course, most such formulas 
focus on the shape between and 1, since only that sec- 
tion is needed for the regular correlation energy. But any 
such approximate formula can be equally applied to K^, 
yielding an approximation for the decorrelation energy. 
Peach et al. [S] analyze various parametrizations for W\, 
and the same forms can be used as well to parametrize 
Kfj,, based on the similar shape of W\ and curves. In 
general, application to will yield a distinct approx- 
imation to the ground-state energy, with quite different 
properties. 
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To give just one example, one of the earliest sensible 
smooth parametrizations is the [1,1] Pade of Ernzerhof 



W x = 



bX 

^cA 



(23) 



One can imagine applying it with inputs of e.g., E x , 
Wq given by Goiiing-Levy perturbation theory, and Woo 
from the SC limit. It yields a sensible approximation to 
W\ in the to 1 range, but because it was not designed 
with the strictly-correlated limit in mind, the formula 
itself is not immediately applicable to the decorrelation 
energy, since, e.g., vanishes. However, much more 

sensible is to make the same approximation directly for 
instead, if one is doing an SC calculation, i.e., 



TABLE I: Comparison of several quantities for three approx- 
imations to Kfj,. Note: ISI uses Ko as an input. The exact 
values are taken from Ref. 151. 





exact 


ISI simp 


exponential 




0.416 


0.416 0.383 


0.456 




0.058 


0.054 0.059 


0.058 


Erie 


0.189 


0.191 0.190 


0.193 



we do an exponential fitting using the method of least 
squares, with the exact values (between /i — 0.5 and 
1) taken from Ref. |5J. We plot vs. /i in Fig. [4] and 
compare several quantities in Table |T] 



K^ = a 



(24) 



whose inputs could be Kq, 



and, e.g., a GGA for 



Kfj, = i. This is then a very different approximation from 
the same idea applied to the usual adiabatic connection 
formula. 

On the other hand, there are several approximations 
designed to span the entire range of A, the most famous 
being ISI (interaction strength interpolation) model [T3] 
developed by Seidl et al. This model uses the values and 
the derivatives of W\ at two limits, namely the high- 
density limit (KS system, A = 0) and the low-density 
limit (strictly-correlated system, A — > oo), to make an 
interpolation. Another approximation to W\ is devel- 
oped in our previous work |25| , which employs Wq , Woo 
and Wq as inputs. We compare the approximate K^s 
obtained from the two models. 

Hooke's atom is a two-electron system (i.e., with 
Coulomb repulsion) in a spherical harmonic well |26| . 
Using the accurate values Wq = —0.515, Woo — 
-0.743, = -0.101 [SJ, and W^ = 0.208 [27], we find: 



K 1 * 1 = -0.947^ + 1.029A/X - 



0.336 



-0.270/ilnB, (25) 



where A = ^/l + 0.653/^ 2 and B = A - 0.263. 
same data substituted in H 7811 " 13 125 , we find: 



With the 



T^simp 



where a = 
Kf, in Fig. 



0.228 



(a 3 - a 2 + 1) + 1.287^(a- 1), (26) 



Jl + 0.354/^ 2 . We plot the two forms of 
4~| The exact curve (down to /i = 0.5) is 



taken from Ref. [8]. We compare three quantities in 



Table 



Although K 1 ^ 1 contains a spurious nln/j, term as 
fj~ |l5lfrfl[23] . it nonetheless yields accurate results. 
The simple model, applied with the usual inputs, is less 
accurate pointwise, but integrates to an accurate value. 

We can try the simple exponential parametrization Eq. 
( 18 1 for again for Hooke's atom. Because we do not 
know the value of d/dfiln at fi = exactly, instead 



0.5 



0.4 



0.3 



0.2 



0.1 



exact 
exact K 
ISI 
simp 
exponential 



0.2 0.4 0.6 0.8 

\1 



FIG. 4: Adiabatic connection curves for Hooke's atom. The 
exact curve (down to /i = 0.5) is taken from Ref. [5J. 



C. H2 Bond Dissociation 

Bond dissociation of the H 2 molecule produces a well- 
known dilemma in computational chemistry 28, 29, 30, 
IBTj . In the exact case, as the bond length R — > 00, the 
hydrogen molecule should dissociate to two free hydrogen 
atoms, with the ground state always a singlet and spin- 
unpolarized. However, spin-restricted, e.g., restricted 
Hartree-Fock or restricted Kohn-Sham DFT, give the 
correct spin multiplicity, i.e. the wavefunction is an 
eigenfunction of S 2 , but produce an overestimated to- 
tal energy, much higher than that of two free hydrogen 
atoms. Spin-unrestricted, e.g., unrestricted Hartree-Fock 
or unrestricted Kohn-Sham DFT, give a fairly good total 
energy, but the wavefunction is spin-contaminated, i.e., 
the deviation of < S 2 > from the exact value is signifi- 
cant. This is known as "symmetry breaking" in H 2 bond 
dissociation. 
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Fuchs et al. [3J] argued that DFT within RPA (ran- 
dom phase approximation) gives a correct picture of 
the H 2 bond dissociation within the spin-restricted KS 
scheme. They also gave highly-accurate adiabatic con- 
nection curves for ground-state H2 at bond length R = 
1.4A and stretched H 2 at bond length R = 5k. The 
curves were interpolated particularly between A = and 
1, shown as the difference of the integrand, AH^, be- 
tween the stretched H2 molecule and two free H atoms 
(Fig. 1 and 3 of Ref. [§2]). 

For R = 1.4A and R = 5k, if we use an interpolation 
(see Ref. 63 in Fuchs paper [35]) to estimate AW\, we 
find reasonable values A Woo = —7.00 and AIUoo = 0.13, 
respectively. Using Eq. (J9|, we find AE DC = 4.96 and 
AE ac = 0.69, respectively. We compare AE DC and AE C 
values in Table [II] The comparison shows a physical ex- 
ample where the strictly-correlated system is a better 
starting point in the calculation. 

TABLE II: Comparison of several quantities for stretched H2 
at different bond lengths. The values for AE X and A_E X c are 
taken from Ref. [32]. All values are in eV. 



bond length 


1.4A 


5A 


00 




-0.98 


5.85 


8.5 


A£xc 


-2.04 


0.82 


0.0 




-7.00 


0.13 


0.0 




-1.06 


-5.03 


-8.5 


AE DC 


4.96 


0.69 


0.0 



We can see that at the equilibrium bond length, |A£? C | 
is much smaller than AE DC , presumably making it eas- 
ier to approximate the ground-state energy starting from 
the KS reference system. This is typical at equilibrium. 
However for stretched bonds, Ai5 DC is much smaller than 
|A_E C |, and so AE DC instead may be better accurately ap- 
proximated in the calculation and the strictly-correlated 
system should be a better reference. Molecules with 
strong static correlation, such as Cr 2 and O3, might fall 
somewhere inbetween. 



IV. CONCLUSION 

In this paper, we constructed an adiabatic connection 
formalism for the strictly-correlated system. We found 
that this adiabatic connection formula and curve can be 
well defined with respect to this new reference. Our for- 
mula connects the strictly-correlated system and the real 
system, by Eq. Q . We also defined the quantity for this 
new integral "decorrelation energy" and related this with 
the usual KS adiabatic connection. We illustrated how 
the decorrelation energy behaves, using the uniform elec- 
tron gas, Hooke's atom, and stretched H 2 as examples. 

We emphasize again that a real application of this 
theory is only possible when the reference, i.e., the 
strictly-correlated system, can be routinely calculated. 
At present, one can calculate quantities such as U sc ex- 
actly only for spherical symmetric systems |14j . However, 
nonempirical approximations to E xc of Kohn-Sham the- 
ory can be employed to estimate Woe with useful accu- 
racy |33j . The computation of this quantity may become 
much easier in the future [34]. If this is so, then based 
on the properties discussed here, the strictly-correlated 
reference may be preferable in cases that are difficult for 
standard KS DFT calculations with standard approxi- 
mations to E xc . In fact, a recent work [34] indepen- 
dently shows progress using exactly the formalism dis- 
cussed here and suggests approximations to E DC . In any 
event, the advent of strictly-correlated calculations opens 
up a whole new alternative approach to DFT calculations 
of electronic structure, and only experience can show how 
and when this will prove more fruitful than the tradi- 
tional (KS) scheme. 
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